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Abstract 

The next-to-leading nonperturbative contributions to heavy quarkonium systems 
are calculated. The applicability of the Voloshin-Leutwyler approach to heavy quarko- 
nia systems for the physical cases of Bottomonium and Charmonium is investigated. 
We study whether the background gluon field correlation time can be considered to be 
infinity or not, by calculating the leading correction to this assumption and checking 
whether the expansion is under control. A phenomenological analysis of our results is 
also performed. The results make us feel optimistic about the T(IS') and to a lesser 
extent about the J/ip but do not about higher levels. We also briefly discuss the 
connection with different models where a finite gluon correlation time is introduced. 
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1 Introduction 



In the early eighties Voloshin and Leutwyler (VL) [111 |2|, |3j developed a theory for quark- 
antiquark bound state systems in the limit itlq — > oo from first principles. In this approach 
it was considered that the quarkonia system can be mainly understood as a Coulomb type 
bound state (this is certainly true if the heavy quark mass is large enough). Moreover, the 
nonperturbative (NP) contributions are corrections and can be included systematically. The 
leading NP corrections were calculated in refs. 0, [|. 

Let us briefly discuss the main features of this approach. First, for calculation conve- 
nience, it is usually chosen the modified Schwinger gauge fixing for the background gluon 
field (see §) 

x l A\x,t) = 0, A°{0,t) = 0, (1.1) 

while the Coulomb gauge is used for the perturbative gluon field. The Hamiltonian then 
reads 



H = H QQ + H g + Hj, (1.2) 

where 

H QQ = P S H S + P 8 H 8 , (1.3) 

A C F a a {fi) A 1 ojjj) 

tis= , #8 = ^7777 > \ LA ) 

m r m Z1S C r 

P s (P 8 ) is the singlet (octet) projector on a quark- ant iquark pair, a s is defined below, H g 

acts on the gluonic and light quark degrees of freedom and 

H I = - 9 -Cx-E a (0), (1.5) 

where £ a = tf — t\ [t\ (t|) is the color SU(3) generator for the quark (antiquark) with 

t® 2 = ^) and E a is the background chromo-electric field. Giving ( |1.5|) for Hj we are 
assuming that the multipole expansion holds. That is 

\x\A Q cd « 1 (1.6) 



but we will also assume Hj to be small as far as we will work within standard quantum 
mechanics perturbation theory (for consistency, we will demand the NP corrections to be 
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small)[]. Thus, ( |1.6| ) becomes 

^»A Q cd, (1.7) 
n 

where j3 n = CpOL s jn. Nevertheless, for obtaining the standard VL results other assump- 
tion is needed, that is, to consider the background gluon field correlation time T g to be 
approximately infinity 

i « 1 — mfa » 1 ~ A QCD . (1.8) 
Now, we have a double expansion 

°((^?) 2 ' ( w )2) - We remark that O 

is more likely 

to fail than (|1.7|) (although, certainly, numerical factors can play a role). The main aim of 
this paper is studying the assumption (|1.8|) and its corrections which, obviously, are going 
to be the dominant NP corrections. 

During several years it was generally believed that Bottomonium and Charmonium sys- 
tems were not heavy enough. Thus, the above mention approach could not be applied. 
Nevertheless, the situation was not so clear for the T(IS') and for the n = 2 fine and hyper- 
fine splitting. Later on, assumption ( |1.8j ) was relaxed and T g was taken into account [|, |6], [7| 
(for a recent study see ||). A detailed study of the T(15') and the n = 2 hyperfine splitting 
was done in the last paper of ref . . Unfortunately, these approaches are model dependent 
since they have to deal with an arbitrary function. Usually, an exponential decaying ansatz 
is used. 

Recently, there have been some attempts to perform a rigorous QCD determination of the 
Bottomonium and Charmonium properties || [K| using only the VL approach. It is claimed 
that a consistent theory for the lower Bottomonium levels (n = 1, 2) and to a lesser extent 
for the J I ip is found. Radiative, NP and eventually relativistic corrections were put together 
for the first time. Some general features for spin independent observables can be inferred 
from this study. The leading perturbative and NP corrections are consistently smaller than 
the Coulomb energy for the T(IS') mass. Likewise, the corrections are under control for 
the T(IS') decay width although the result is quite dependent on the scale. For n = 2 
observables the situation is much more doubtful since the perturbative and NP corrections 
are almost as important as the leading term. Nevertheless, quite good agreement with the 
data was obtained. For these observables and the T(IS') decay width the NP corrections 
play a fundamental role in order to agree with experiment. 



x If we did not consider Hi nt to be small, \x\ would be an arbitrary function of toq, a s and Aqcd- This 
more general situation is beyond the scope of this paper. 



4 



It could be argued that further orders in the relativistic, radiative and NP corrections are 
needed in order to improve these results. Moreover, so far, no model independent check of 
approximation ( [L.8|) has been done though it is the most likely to fail. This check is clearly 
very important to clarify in which situation we are, that is, whether ( |l.<j| ) holds or not and in 
which observables. In the former case the VL approach should be applied while in the last 
case T g has to be taken into account in some way and the models displayed in refs. f|-|| 
may be applied. In order to disentangle this issue, we calculate the next-to-leading spin 
independent NP contributions to the energy levels and wavefunctions of a heavy QQ pair. 
We obtain analytical expressions for them, and we check whether the expansion is under 
control. 

Let us briefly comment upon the fine and hyperfine splitting (n = 1,2). For them the 
agreement with the data (when available) is quite good according to ref. |H| (although errors 
are large). This is quite surprising since Dosch and collaborators J7| also find agreement with 
the data using a model where T g is taken into account. Undoubtedly, this should require a 
detailed study which goes beyond the scope of this work and we will not perform it here. We 
expect these observables to depend on shorter distances being the NP contributions smaller. 
Combinations of these observables can be built such that the NP corrections are under 
control as it was done in JTTJ (see also |12| where a omitted contribution was calculated). 

We distribute the paper as follows. In sec. 2 we calculate the next-to-leading NP correc- 
tions to the energy and decay width. In sec. 3 we perform a phenomenological study of our 
results. The last section is devoted to the conclusions and the discussion of models where T g 
is taken into account. A few formulas are relegated to Appendix A and we discuss the role 
played by the octet potential in Appendix B. 



2 Theoretical results 



We display the calculations and results for the next-to-leading NP corrections to the spin 
independent energy levels and wavefunctions. 



2.1 Energy Levels 

In this subsection we calculate the next-to-leading NP contributions to the energy levels. 
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Using the multipole expansion (17?) and standard Quantum Mechanics perturbation theory 
the energy correction reads 



5E nl = {0\(n,l\Hj- ^ Hl \n,l)\0) 

= ^(0\Ef(0)(nJ\r- ± -r\n, l)E}(0)\0) . (2.1) 

18 E n - H 8 - tig 

E n is the Coulomb singlet bound state energy, n, I are the radial and orbital quantum 
numbers. As far as we do not study the fine and hyperfine splittings the corrections do not 
depend on j (total angular momentum) and s (spin) so we will not display these indices in 
the states. H g ~ Aqcd and H$ ~ E n . 

The octet propagator mixes low 0(Aq C d) and high energies 0(E n ). No assumption 
about (|1.8| ) has still been done in (|2.1|) . Let us assume ( |1.8|) holds, it follows that an OPE 
can be performed where the parameter expansion is of order 



H g \ 2 ( Aq C d\ 



2.2) 



and we obtain 



where 



E n -H 8 J \mPl) 

oo 

6E nl = J2C r O r , (2.3) 



r=0 



/ 1 \ 2r+1 

C r = M ff { En _ H J r >'0 ( 2 - 4 ) 

O r = ^(0\E«H?E«\0). (2.5) 



and 



18 

Using the equation of motion, the gauge fixing and Lorentz covariance we obtain 



a 2 

O r = -—v Po ...v^v ao ...v ar (2.6) 
54 

x(0|Tr([D ft (0),[...[Z^(0),^ 



where v is the velocity of the center of masses frame with v 2 = 1 (in the comoving frame 
v = (1, 0)) and the trace is in the adjoint representation. 
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Let us give the order of each term (< r >~ " 



m(3„ ■ 



< r ) = C,0, ~ mQl ( ^^-) 



2 / . x 2r+2 

l^QCD 



For 5E^i we obtain the standard VL result [|[ [|. Nevertheless, we are interested in the 
next correction 

= <7 l0l ~ m#%|£- , (2.8) 



where 

Ci = m? (e^wS (2 ' 9) 

After a rather lengthy calculation we obtain 



<7 1 = — ^lf(n,0, (2-10) 

where if (n, /) is a dimensionless function which can be obtained using techniques developed 
in ref. @. We refer the reader to Appendix A for the analytical expression which can be 
written as a rational function. There, we also provide some numbers for the lower levels. It 
can be easily checked that it has the right leading n 10 dependence. 

For Ox, using Lorentz covariance, Bianchi identities and the equation of motion, we 
obtain 

Oi = I -S 4 £<°l = q A t a W A q B t a rq B : |0> + ^<G 3 > \ , (2.11) 



(G 3 > = g 3 f abc {0\ : GW : |0> . (2.12) 



where 

/r<3\ — „3f _ /ni . r 

and A, B are SU(3) flavor indices. 

We can work (|2.11|) further using SU(3) flavor symmetry and the factorization hypothesis. 
Finally, it reads 

0i = ^^WMmW + \^))- (2-13) 

The physical cases will be studied in the next section. Anyway, let us comment that 0\ 
is going to be positive and, therefore, the energy correction to be negative. 

In Appendix B the error made if we neglect the octet potential in (p.9|) is discussed. 



7 



2.2 Decay Width 

In this subsection the corrections to the decay width of the n 3 Si levels of quarkonium 
are calculated. For further details we refer to refs. |l|, 0. 
The decay width of the n 3 Si levels of quarkonium reads 



T(n 3 S 1 -> e + e-) = it 



_M(n 3 Sx) 
where Q is the quark charge and 



... , , , (x = 0\G s (E)\y = 0} 1 (2.14) 

07T / "'-"'pole 



G S (E) = P s (0\j^\0)P s (2.15) 

is the nonrelativistic propagator (Green function) of the system projected on the colorless 
sector of a quark-antiquark pair and the gluonic vacuum. E is the energy measured from 
the threshold 2m. Near the pole n we have the expansion 

(x = 0\G s (E)\y = 0) = - P 2^ Pn + Q((E n + 6E n0 - E)°) 

h, n + dh, n0 — hi 

= ~ { e!-EY + E~^E + ° {{En ~ E)0) + ° ( ^ o) • (2 ' 16) 

On the other hand using ( |1.7|) we obtain 

(x = 0\G s (E)\y = 0)~(x = 0\G^(E)\y = 0) + (x = 0\5G s (E)\y = 0) , (2.17) 

where 

(x = 0\Gf\E)\y = 0) = (x = 0\-^—\y = 0) = — + 0((E n - E)°) (2.18) 

ri s — hi hi n — h, 

_ 1 mPn x 3 

71 I 

and 

(x = 0\5G s (E)\y = 0) = P s (0\(x = 0|^4-^y 7 — -L__^_l_|^ = )|0)P S 

/V ^'"" • t^— • 0((J$, • (2.19) 



(E n — E) 2 E n — E 
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Proceeding in the same way than in the preceding section we are able to split low from high 
energies (OPE). 

oo 

(x = 0\5G S {E) |j? = 0) = £ °rO r , (2.20) 

r=0 

where 

1 / 1 \ 2r+1 1 

C? =(f=o| f f \y = o) 

V l H 8 -E \H 8 -EJ H s -E iy 1 

Ar) Ar) 

' 2 + j- + Q((E n - E)°) (2.21) 



(E n - Ef (E n - E) 

and O r is the same than above. Now, from these expressions we can read off the observables 
we are interested in, namely 



*Pn = £ ^S = £ 430r , (2.22) 



r=0 r=0 
.1 00 



6E„ a = — Y, A -lOr. (2.23) 

Pn r=0 

This provides a new method to obtain the energy corrections for I = states. We have 
used them in order to check our results of the previous subsection. 5p^ and 5E^q were 
already calculated in ref. 0. We have used their results to check ours for this lower order. 
We are mainly interested in the next correction Sp^ which, after some effort, has been 
calculated exactly. We write it within the dimensionless quantity 

where the analytical expression and numerical values for W(n) are displayed in Appendix 
A. We can see that it has the expected leading n 12 dependence. 

In Appendix B we give analytical expressions for 5p^ and SE^f (I = 0) with an arbitrary 
octet potential, | — ► Ar2 1 _ 1 . There, we study the error made by neglecting the octet potential 
(N c — > oo) in our results. 



3 Phenomenological analysis 



In this section we confront our results with the data. We expect our corrections to be 
small when the theory works and the opposite could be a signal that the whole approach 



9 



breaks down. We will try to see the effect of the new contributions on the results found by 



the authors of ref. M, 10 1 and we will mainly follow their notation. 



Let us briefly discuss the set of parameters we are mainly going to use here. We make 
Aqcd running between a range of values compatible with those used by the authors of ref. 

A n Q ^ = 250tT Mev. (3.1) 

The mean value approximately corresponds to the value found by Voloshin in ref. using 
sum rules. The upper bound corresponds to the standard value of Aqcd found in DIS (Deep 
Inelastic Scattering). The lower bound is closer to the values proposed by sum rules some 
years ago. For the two gluon condensate we quote the value given in ref. [f§ 

(aG 2 ) ~ 0.042GW 4 . (3.2) 

We will allow 0\ to have some errors. It is not our aim to give rigorous errors but to see 
whether the results are sensitive to sensible variations of 0\. For the three gluon condensate 
we choose the value 

(G 3 ) ~ 0.045 ± 0.009GW 6 , (3.3) 

which comes from the dilute instanton gas approximation [|TJJ] while the errors have been 
estimated using the value obtained in ref. fT5| . Similar results are obtained using the three- 
gluon condensate values given in ref. |16| (sum rules) or ref. |7| (lattice)^ 
In the sum rule approach it is usually found the quantity 

K = a s < 0\qq\0 > 2 (//) , (3.4) 

which appears to be very weak scale dependent once the anomalous dimension is taken into 
account so it is usually considered to be a constant. In principle this expression could be 
calculated once the quark condensate and a s are known at some scale. However, there 
has been some controversy about the error made within the vacuum saturation hypothesis 
and that may be the four-quark condensate has been underestimated by around a factor 2 
18l-||20||. We just take the average value of these references given in ref. |21 



k = 3.8 ± 2.0 x 10 GeV . (3.5) 



2 Notice that < G 3 > calculated in euclidean space has an opposite sign compared to that in Minkowski 
space. 
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Nevertheless, instead of Eq. ( |3.4| ) we actually have 

a] < 0|gg|0 > 2 (ji) = a s {fi)K (3.6) 

in 0\. The only remaining scale dependence is in a s {fi). We stress that a s has to be computed 
at the subtraction point scale, generally, the inverse Bohr radius. This value depends on the 
physical system one is studying. Therefore, 0\ is going to be scale dependent, although 
weakly. So finally, we obtain using (|3.3|) , Q3.5[) and their errors (slightly larger errors would 
not change the conclusions) 

O x {n) = (3.13±g;|| + a s (fi)7AltH) x lO^GeV 6 . (3.7) 



Let us finally remark that the range of values given in Eq. ( |3.1| ) does not include those 



obtained from LEP and r-decay date pl[-||23||. We quote the value from reference ref. p4 

a s (M z ) =0.118 ±0.003 



which approximately corresponds to 



^QCD — 



383 



1-52 
-48 



(3.J 



(3.9) 



There has also been some controversy about the value of the gluon condensate < a s G 2 > 
and several claims of larger values can be found in the literature. See [24] and references 
therein where the average value 



< a s G >- 



0.081Ge1/ 4 



(3.10) 



was given excluding the SVZ-like value [fL4]| . In fact, for many of these studies the input 
values of Aqcd are compatible with those given in Eq. ( |3.9[ ). This region of parameters was 
uncovered in ref. || |10[ . We have also studied these new possibilities although our principal 
aim will be to compare with the results of ref. 0, |1O[0. I remark, in order to be consistent, 
that < G 3 > should be also changed since |T7]] 



< G 3 > oc < a q G > 



(3.11) 



3 A complete analysis of this new range of parameters, including the static two loop potential and the 
relativistic corrections, is under way [p5| . 



11 



In our case (changing accordingly the errors) 

< G 3 > — >< G 3 >= 0.087t° ° llGeV 6 
and (using the value of k given in (|3.5|)) 

1 ( f x) — > OM = (6.041^ + a s (fJ,)7AltH) x lO~ 4 GeV 6 



(3.12) 



(3.13) 



Now we have two possible sets of parameters. We choose Eqs. ( |3.1|) , ( |3.2|) and (|3.7f) to be 
our first set of parameters or set I, while Eqs. ( |3.9| ), ( |3.10| ) and ( |3.13| ) are our second set of 
parameters or set II. 

Let us now write the general mass formula for any state 



where 



M{n, I) = 2m b + A 2 {n) + A 3 (n, I) + SE, 



(0) 
nl 



5E. 



(i) 

nl 5 



(3.14) 



A 3 {n,l) 



-2m, 



A 2 {n) = - 

C 2 f p a 2 s (n)a s (fi) 
8irn 2 

6E$ = m b 



2m b 
( In 



in" 



fin 



rn b Cfa s (n) 
e nl n 6 ir(a s G 2 ) 



+ ip{n + l + l) 



(3.15) 



(m b C f a s ([i)) 4 ' 

e n i were first calculated by Leutwyler in ref . . For the levels we are interested in they read 



1872 



2102 



9929 



£io 



^20 



£21 



(3.16) 



1275 ' "~ 1326 ' 9945 
a s is the two- loop running coupling constant. A 3 , which include radiative corrections, and 



l + (ai+7s/3o/2)- 



7T 



(3.17) 



a i 



31C A - 20T F n f 
36 



have been quoted from ref. 0. We remark that the relativistic corrections have not been 

included in this analysis. This is due to the incomplete knowledge of the 0(ma 4 ) corrections 

for the mass (or the equivalent for the decay width) in the MS scheme because the static two 

loop potential has not been calculated^. Although naively, these corrections are expected 

4 Recently, we became aware that the static two loop potential had been finally calculated ^6|. Wc 
expect to introduce those results in the near future p5fl. 
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to be small they could blow up if a s , a s are big enough. In order to make quantitative 
this statement the next-to-next-to-leading perturbative contribution, including radiative and 
relativistic corrections, should be calculated. Moreover, it has been seen, according to the 
recent sum rule derivation of the Balmer mass formula [ 23 1 , that the relativistic corrections 
tends to compensate the radiative coulombic corrections. 

Let us briefly comment upon some general features of the numerical study. Large values 
of a s , a s appear in the problem. Moreover, for larger values of Aqcd the value of a s also 
increases. Therefore, for the second set of parameters, the results are going to be more 
doubtful. This is especially so for n = 2 Bottomonium states and Charmonium. Neverthe- 
less, without the knowledge of the next-to-next-to-leading perturbative contribution, we do 
not have a real check in order to know when the a s perturbative expansion breaks down. 
Therefore, we refrain from displaying numbers for n = 2 Bottomonium states and Charmo- 
nium with the second set of parameters here (although we do not rule out the possibility the 
VL approach to work for them). This is relegated to future work where the static two loop 



potential and the relativistic corrections will be taken into account 25 



3.1 T(15) mass 

We start our study with the T(1S) mass where the formalism is expected to apply better. 
We proceed as follows. First of all, we fix mj and the inverse Bohr radius a bb \ from the self- 
consistency equation 

-i rn h Cfa s {a^ x ) 
a bb,i = 2 ( ) 

and the T(1S) mass set at the inverse Bohr radius scale (/i — > %, x ). We allow for the both 
above mentioned sets of parameters (without using errors for 0±) and give the relative weight 
of each contribution in table I. It can be seen that Aqcd turns out to be the main source of 
error. We also get that perturbation theory works nicely here. Similar results than in ref. 
|| are found being not to much affected by the new contribution. 

For the first set of parameters, plugging the errors given in Eq. ( |3.7|) into dEy^ and 
making the ratio with 5E$ , we find 



• ,5£fo\ 
0.19 < — ^ <0.38, 0.14 < — ^ <0.27, 0.10 < — £j < 0.20 , 3.19 

SE® 6E® 6E® 



for Aqcd = 200, 250, 300 MeV respectively. Likewise, for the second set of parameters, we 
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A n/=3 

IV QCD 


A 2 


^3 






m fe 


a bb,l 


250H8 


-376$ 




is;? 


-4±| 


4881±|f 


1355l 12 i 


QOQ+52 
000_ 48 


doo +68 


93±i| 


17+5 


-11? 


4943^ 


1626±iif 



Table 1: We display A 2 (l, 0), A 3 (l, 0), c^i? and 5E$ for the T(15). The last two columns 
give our results for m b and a^ b \. All the quantities are displayed in units of MeV. We have 
used f^QCD = 330l|g MeV for the second set of parameters. 



obtain 

5E {1) 5E {1) 5E {1) 

0.08 < I — ^|<0.14, 0.06 < I — £jcl<0.11, 0.05 < I — ^|< 0.09, (3.20) 

6E$ 6E[$ 5E$! 

for Aq-^ = 282, 330, 382 MeV respectively. We see 5E[^ is always reasonably smaller than 
5Eiq and we can safely conclude we are inside the VL regime for this observable. 

Let us comment on some features of the results. They are very stable under variations 
of \i because the \i scale dependence of A 2 and A 3 cancel each other. This can be seen in 
more detail in Figure I where we plot the T(1S) mass as a function of the scale. Moreover, 
the scale of minimum sensitivity and the inverse Bohr radius scale are quite close. Indeed, 
the difference between the Mmsn values for those two scales is completely negligible. Other 
very important point is that not only the NP corrections are small but also under control. 
That is, higher order NP corrections can be calculated giving smaller contributions. The 



values we obtain for mj are compatible with the lattice calculation 27 although somewhat 



larger than those obtained in sum rules |28 if the comparison is done with the two-loop 



perturbative pole mass. Nevertheless, the agreement is quite good if the comparison is done 
with the three-loop pole mass. 

Let us finally point out one prediction which follows directly from our results, A, the 
nonperturbative parameter relating the mass of the heavy meson (B, D) to vhq (nib, Tn c ) in 
the HQET We obtain 

A = 433+%MeV {A n Q f cn = 260™ MeV) ; A = 370-f 8 MeV (A^= D 4 = 330^MeV) , 

(3.21) 

for the first and the second set of parameters respectively. Similar comments than for m;, 
apply here when comparing with the literature. 
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FIGURE 1 



9.495 



9.485 



9.475 



9.465 



9.455 




Figure 1: Plot of My(is) (GeV) versus \i (GeV). The NP corrections have a little effect on the 

draw. We have used the first set of parameters with J^cd = 250MeV r and = A881MeV. 

V= 3 



The same features hold for the other values of A 



qcd- 



3.2 n = 2 states 

We now discuss n = 2 states and their masses M(2, 0) and M (2, 1). We will only use the 
first set of parameters in this subsection. We expect the theory to work worse or even to fail 
completely here. In fact, there are several signals that for n = 2 the theory does not really 
work. For instance, we can now predict the physical masses since we have got rrib and, in 
order to set the scale, we can use the self-consistency equation for n = 2, that is 



m b C f a s (a b ^ 2 ) 



l bb,2 



(3.22) 



nl I 



> \A 2 \ and also \8E^ 



nl 



> 



We display the results in table II and III. We find out that \SE, 
\8E^i\, excepting the case| n = 2, 1 — 1 for ^qcd = 300MeV A , so the theory is not really 
trustworthy. Agreement with the data is also very bad. We have also studied the dependence 



Nevertheless the reliability of perturbation theory is marginal since \8E\[>\ < \8E^'\ and \5E^>\ > ||A 2 |. 
Moreover, a s — 0.403 and a s — 0.628 are quite big making the a s expansion rather doubtful. Therefore, 
knowledge of perturbative and NP next order contributions would be welcome in order to discern whether 
the a s and/or the NP expansion work or not. 
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ly QCD 


A 2 


^3 


°- c '20 


5E% 


M(2,0) 


a bb,2 


250118 


-178+1 


-29;? 


ooq-109 






qoo+095 



Table 2: We display A 2 (2, 0), A 3 (2, 0), and JE^. The last two columns give our results 
for M (2, 0) and a^ 2 . All the quantities are displayed in MeV. The experimental value is 
M(2,0) = 10023Me1/. 



iv QOD 


A 2 


^3 




6E$> 


M(2,l) 


25011° 


-i78;l 


_ ?1 -19 


91 0-O68 
z±o +120 


-2471™ 


94791^ 



Table 3: We display A 2 (2, 1), A 3 (2, 1), (5^ j and The last column gives our results 

for M(2, 1). All the quantities are displayed in MeV. The experimental value, averaged over 
angular momentum, is M(2, 1) = 9900MeV. 

on the scale. We find that the scale we reach the minimum and the scale we obtain doing 
the self consistency equation (with n = 2) are not so near now, especially for the n = 2, 
I = 1 state. Therefore, we rather try with mass shifts like 

AM(2,0) =M(2,0) -M(1,0), AM(2, 1) = M(2, 1) - M(l, 0) , (3.23) 

as it was indeed done in ref. ||. Let us briefly comment upon this work. The results found 
there were already quite delicate because the corrections were around to be as large as the 
leading term. However, this fact was compensated with the good agreement with the data 
obtained. 

AM(2, 0) = 479MeV (expt. 563MeV) , AM (2, 1) = All MeV (expt. 450MeV) . (3.24) 

From this point of view we should say that the new contribution completely destroys 
this agreement as we can see if we simply plug our result with their input set of parameters 
(m b = 4906MeV, \i = 986MeV for AM(2,0) and /i = 1062Me1/ for AM(2, 1)) and with 
the errors given by ( |3.7|) . 

AM (2, 0) = -182l|£MeV , AM (2, 1) = -21V^MeV . (3.25) 

The new contribution even changes the sign of AM! This result is due to l&E^ 1 )] > l&E^I. 
It can also be seen that the dependence on the scale is very strong. Let us briefly explain 
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what is going on here. Let us first neglect 5E^ . Then, for consistency of the theory, one 
demands 

L4 2 |>L4 3 | \M > \6E®\. (3.26) 

The first constraint does not allow to lower the scale while the second one does not allow to 
increase it, hence, only a very small window is permitted. fl3.24|) was found for a /x scale such 



as radiative and NP contributions were equal in absolute value and even in this optimum 
situation ( p.26|) was not well satisfied (\A 2 \ ~ \A 3 \ ~ \SE^\). Now, if we introduce the 



new contribution a new consistency constraint |5i?( )| > arises which is also not well 

satisfied or, even more frequently, it is not satisfied at all. 

We do not give more numbers here as far as we conclude the theory is not really trustwor- 
thy for n = 2 states. Indeed, all these results make us feel that probably both approximations 
( [L.7D and ( |1.8|) fail. In fact, that can be seen by the large mass gap between the n = 1 and 
n = 2 states which can not be obtained by only perturbing about the Coulomb spectrum. 
Therefore, n = 2 states seem to behave completely different than n = 1 states. For the for- 
mer the NP contributions are large and can not be treated as corrections so other approaches 
should be attempted. For example, a hopeful approach has been developed by Dosch and 
Simonov [|30| where they seem to be able to connect the VL regime with the regime where 



potential models work. 



3.3 Charmonium 

We now briefly discuss the Charmonium case. We will only use the first set of parameters. 
For this particle the Bohr scale is rather small and consequently a s quite large (indeed larger 
than for n = 2 Bottomonium states) so the a s perturbative expansion should be taken with 
care. Let us anyway say a few words about it. First let us find the value of our contribution 
with the set of parameters given in ref. || (m c = 1570MeV and /i = lGeV). We obtain 
(using the errors given in ( |3.7p ) 

5E[l } = -154+llMeV, (3.27) 

\5Eiq | < |5-E'io' , | and \A>z\ — ISE^l (this was the condition used to fix the scale in ref. ||) 
for those parameters so the results are a little bit doubtful. Let us now use our standard 
procedure by fitting the scale with the self-consistency Bohr radius equation (notice that we 
use 7] c rather than J/ip in our procedure). We display the results in table IV. Everything 
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ly QCD 


A 2 


^3 








a cc]l 


250118 


-368+so 


781^ 


KE-1B 
dd +38 


-1411 


161411 


770lp 



Table 4: We display A 2 (l, 0), v4 3 (l, 0), 5E^ and <5-E , io ; for the n c . The last two columns give 
our results for m c and a~\. All the quantities are displayed in MeV. 



(i) 



seems to behave quite properly, especially the NP contributions, but for those scales a s is 
quite big so it would be desirable to obtain the next-to-next-to-leading order correction 
in order to discern whether the a s expansion expansion can be applied or not. We notice 



here that 8E^ and SE^ behave quite well now, say, the inequalities |v4 2 | > l^-Eiol > l^io 



If I > \se$?\ 

1 and the blowing up of the NP corrections 



are well satisfied. The reason being that n 
with n does not arise. 

Just as in the m;, case the values we obtain for m c are somewhat larger than those 



obtained in sum rules [28[ if the comparison is done with the two-loop perturbative pole 
mass. Nevertheless, the agreement is quite good if the comparison is done with the three- 
loop pole mass. 



3.4 T(15) decay width 

Let us conclude with the T(15) decay width. 



r(T(l£)) = r<°> (1 + (Sr) (l + <W + ^ + ^-l , (3.28) 

2pi 2pi 



where 

r (o) 



Q b a 



1 2 



em 



m b C f a s (p)) , 6 WF = 3f3 In — - 7£ 



M T (is)\ V \rn b Cfa s (fi) J J An 

AC f a s (fi) 5pf ] _ / 270459 1838781 \ tt < a s G 2 > 
r ~ tt ' ~2p7 ~ V217600 + 5780000/ mfa%fi) ^' ' 

We have quoted ( |3.28|) from ref. || adding our result. Let us first discuss the set I of 
physical inputs. We have drawn F with and without Sp^ in Figure 2. For small values of p 
the results are not reliable because the perturbative corrections become too large. It can also 
be seen the strong dependence on the scale. The results on Figure 2 should be compared 



FIGURE 2 
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1.5 




0.0 1 1 

1.4 1.8 2.2 2.6 

n 

Figure 2: Plot of the T(1S) decay width (KeV) versus (i (GeV). We have used the first set 
of parameters with F^qcd = 250MeV r and mj = 4881MeV r . The continuous line draws the 
decay width with 0\ = OGeV 6 , the dot-dashed line with the mean value displayed in (|3.7|). 
the dashed and dotted line, respectively with the lower and upper value displayed in (|3.7|). 
The same features hold for the other values of ^qcd- We also plot the experimental value 
r(T(LS)) = l.UKeV. 

with the experimental value r(T(15)) = lMKeV and with r(T(15)) = l.l2KeV ', the 
result found in ref. |TU] using a value of \i ~ 2.33GeV A {nj = 4) such that the perturbative 
and NP corrections cancel each other. We now include our result for such a scale. We find 
(using the errors given in ( |3.71 )) 

r(T(l£)) = UMt^KeV . (3.30) 

We see that it destroys agreement with the data. What is going on here is that for such 
a scale ^ 1 and \8p^\ > \Sp^\ so the result is not reliable. Indeed, one can not get 

agreement with the data for any scale as it can be seen in Figure 2. If one lowers the scale 
in order to fulfill > |<5pi^| the perturbative contributions grow in size with negative 

sign and make agreement with the data worse. 
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These results are somehow surprising since we would expect to obtain the same features 
as those obtained in the Myns) case. One point is that here we are faced with a two scale 
problem, the annihilation scale 0{m b ) and the wavefunction scale 0(a bb x )0. For this kind of 
observables factorization holds, that is, we can split the observable in two pieces the behavior 
of which is dictated by long and short distances respectively. Effective theories are especially 
useful in these situations. In our case the suitable effective theory is NRQCD |3l| which has 
proved to be extremely successful in separating low from high energies. 

Using NRQCD it can be seen that the right scale for 5 r is p ~ mi,. The NP contributions 
behave properly for their natural scales p ~ that is, the inequalities \pi\ > \Sp^\ > 
\Sp^\ are well satisfied, especially the first one. The typical scale of and 5wf is also 
p ~ % bl . The problem is that Swf is large and strongly dependent on the scale. For instance, 
we obtain (with Aj^ = 250MeV, m b = A881MeV and a bb \ = 1355MeV) r(T(15)) = 
Q.UKeV for T^(a bb \) = 2.90KeV, 



r (0) r (1) 

<WKw) = -0-81 , 6 r (m b ) = -0.31 , -^{a bb \) = 0.07 , -^-(a^) = -0.02 . 

ipx zpi 



while for 5wF(2a bb \) = —0.08 and the remaining set of corrections kept fixed we obtain 
r(T(15')) = 1.84KeV. The experimental value lies in between for the value 5wF(2.21GeV) = 
—0.23. Although 5wf is rather large it is not senseless to suppose that higher orders could 
improve the result changing significantly the value of 8wf- 

For the second set of parameters the same general features are observed. The NP con- 
tributions are under control ((|pi| > I^Pi^l > I^Pi 1 "*!) but the a s perturbative expansion 
behavior is even worse. For instance, we obtain (with Aj^ = 330MeV, m b = 4943MeV 
and a bb \ = 1626MeV) T(T(15)) = 0.09KeV for (a^) = 5.02KeV, 

SwF(a bb \) = -0.87, 5 r (m b ) = -0.36, -^-(a^) = 0.05, -^(a^) = -0.01. 

zpi zpi 

Summarizing we may conclude that the r(T(lS')) belongs to the VL regime since for scales 

p ~ a bb \ the NP contributions are under control (\p\\ > \Spi\ > \8pi\). However, no 

precise determination of r(T(lS')) can be given due to the strong dependence on the scale 

and the large perturbative corrections involved. The solution could be to calculate the 

two loop Q — Q potential and introduce the next order perturbative correction. This may 

significantly improve the value of Swf- 

6 The Mx(is) does not have this problem because it is an asymptotic feature of the physical state and 
only the typical bound state scale, a bb j, can appear in the dynamics. 
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From the whole set of results we can conclude that the T(1S) belongs to the VL regime. 
Therefore, T g is not needed and the approach developed in ref. [[7) or potential models should 
not be applied here. 



4 Discussion and Conclusions 

Let us briefly comment upon models where T g is fully taken into account ]7| |8|. There, 
it is obtained 

7r < a s G 2 > , , 1 i 7 \ i a -, \ 

ok/r,! — (n, I r, ; — rAn, I). 4.1 

nl 18 V ' 1 H 8 - E n + 1/T g 41 ' 7 V ; 

In order to compare it with our results we should take the limit T g — > oo in ( |4.1| ). That is, 
1/T 9 << m(3 2 . We see that it has the right asymptotic limit, but it fails to describe properly 
the preasymptotic corrections. For instance, instead of 8E^ ~ ( [4. 1|) already has 



Aq CD corrections. This is not so strange if we recall the exponential approximation for the 



gluon correlator (Euclidean space) which has been used to obtain ([0 



< G(t)G{0) >~ e- t/T v (4.2) 

is only expected to provide a good approximation for long distances which is indeed the 
opposite limit that we are studying. That is, we may only expect ( |4.1| ) to provide a good 
description of the leading NP corrections when ( |1.7|) is satisfied and[] 



T 



» 1710, 



2 



Indeed, lattice calculations suggest that the right behavior at long distances is 



Therefore, we conclude that (pi~T| ) can not be applied when the approximation ( |1.8| ) is satisfied 



which is the situation we are interested in. 
Let us now summarize our conclusions. 

We have calculated the next-to-leading NP corrections to the energy and wavefunctions 
for heavy quarkonium systems with general quantum numbers n and I. We have obtained 
general analytical expressions. We have also been able to obtain analytical expressions for an 
arbitrary octet potential. We have also exactly calculated the error done if we had neglected 
the octet potential (N c — > oo) in our results. We have seen these observables to depend on 
the three-gluon condensate. 



''in this situation new NP contributions seem to appear |32 
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We have reanalyzed the applicability of the VL approach for the lower levels of Bot- 
tomonium and Charmonium. We have studied two possible sets of parameters, especially 
the one compatible with the values used in || |10| . We have carefully studied the multipole 
expansion ( |1.7[ ) and the approximation ( |1.8| ), especially the last one. As we have remarked 
throughout the paper, this last point is essential in order to distinguish whether a finite T g 
is needed. We also believe our work clarifies whether the a s perturbative expansion can be 
applied. 

It has been proven that Mx(is) lies inside the VL regime since the NP corrections are 
under control and very small. It follows that neither Q4.ll) nor potentials models should be 
applied here. Therefore, this observable could provide one of the best model independent 
determination of the rrib mass where the principal source of error will come from Aqcd 
and higher orders in the perturbative expansion which are calculable and do not introduce 
any new parameter. Nevertheless, these corrections could be large. Other very striking 
consequence is that Mx(is) should not be used fixing parameters in the potentials models. 
This would obviously change the value of these parameters. It would be rather interesting to 
see how important these changes are. Finally, we have also seen the weak /i scale dependence 
of the result. 

For the n = 1 Charmonium level, using the first set of parameters, the theory seems to 
behave properly and the VL approach to work but the largeness of a s makes the results less 
reliable. 

For the first set of parameters we have also shown that the VL approach is not really 
applicable for n = 2 states and probably both assumptions the multipole expansion ( |1.7| ) 
and the approximation ( |1.8|) fail for those states and obviously for higher levels. Therefore, 
in order to control the NP contributions, we should look for other approaches, perhaps in 
the spirit of ref. ||30|| , where the NP contributions are not considered to be corrections. 

The T(IS') decay width is very interesting. We have seen that calculating naively we 
do not obtain agreement with the data. The reason is that the new NP contribution blows 
up before we can reach the scale used in ref. Jl0[. It has been noted that the T(IS') decay 
width is a two-scale problem and, in the spirit of NRQCD, different scales should be used for 
the different contributions. We have argued about whether agreement with the data can be 
obtained. Nevertheless, errors are large owing to the large value of Swf for scales /i ~ . 
The solution could be to calculate the two loop Q — Q potential which could fix the value of 
Swf- This may significantly improve the agreement with the experiment. We have also seen 
that the multipole expansion (|L7|) and the approximation (|L8| ) applies for scales [i ~ 
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((|TT7|) is especially well satisfied). 

We have not performed the calculation for the fine and hyperfine splitting but we believe 
that they will behave in the same way than the T(IS') decay width, being strong scale 
dependent. 

Improvements of our results would come from diminishing errors in the parameters. It 
would also be very interesting to have the two loop Q — Q potential in order to obtain 
the next-to-next-to-leading perturbative contribution to masses and decay widths. It would 
improve the rrib results and the T(1S) decay width. As we said, for the latter it could be 
essential in order to set the right value of 5wf- For the r] c and n = 2 Bottomonium states the 
knowledge of the next-to-next-to-leading perturbative contribution could eventually discern 
whether the a s expansion works or not. Moreover, a complete analysis with the second set of 
parameters remains to be done for these states (the conclusions obtained above could change 
for this set of parameters). This analysis will be only reliable once the above perturbative 
contributions are calculated. 

When available, our results may eventually fix parameters in NRQCD and HQET. For 
both of them we can fix rrib and m c in a model independent way. Once the masses are known 
we can obtain the HQET parameter A as it has already been done above. We also expect 
to be able to give QCD rigorous predictions of NRQCD matrix elements related with the 
T(IS') decay width. Moreover, this theoretical framework is able to deal with octet matrix 
elements in a natural way so further predictions may, in principle, be given. 
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A Formulas of sec. 2 

In this section we display analytical expressions and numerical values for the formulas 
obtained in sect. 2. 



1024 n 8 

H (n, I) 



81 (9 n 2 - 64) (81 n 2 - U) 6 (81 n 2 - 256)' 
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x (^13295844358881280 - 68153404341354496 n 2 + 145600287615221760 n 4 

- 168309372752363520 n 6 + 114216987240880128 n 8 - 46158344158975776 n 10 
+ 10789755579716526 n 12 - 1327743092409993 n 14 + 65241222927111 n 16 + 
I (32604849090592768- 128931229385883648 n 2 + 198763892651851776 n 4 
-154057684466810880 n 6 + 65418421737648384 n 8 - 15363511613472780 n w 
+ 1857043938050595 n 12 - 88965303991515 n 14 ) + 

I 2 (40459485281517568 - 155542385042915328 n 2 + 231200955030306816 n 4 
-171789196119588864 n 6 + 70180300981986048 n 8 - 15974396240527980 n 10 
+ 1886699039381100 n 12 - 88965303991515 n 14 ) + 

Z 3 (15799724393103360 - 53517313156055040 n 2 + 65200816249896960 n 4 
-35606723591467008 n 6 + 9549421221866688 n 8 - 1223181401792805 n 10 
+59310202661010 n 12 ) + (A.l) 

1 4 (8125992224686080 - 27496161183006720 n 2 + 33417136857415680 n 4 
-18162612510511104n 6 + 4838867443911744n 8 - 615121070102415 n 10 
+29655101330505 n 12 ) + 

1 5 (271356033761280 - 885005525975040 n 2 + 980074478960640 n 4 
-431 100857733120 n 6 + 76988199574080 n 8 - 4236443047215 n 10 ) + 
Z 6 (90452011253760 - 295001841991680 n 2 + 326691492986880 n 4 

-143700285911040 n 6 + 25662733191360 n 8 - 1412147682405 n 10 ' 



„, , 141912051712 . , 484859657191424 /A , 

if (1 = , H(2, = , A.2 

v ' ; 844421875 ' v ' ; 2040039729 ' v ; 

_ 102150951135870976 _ 1299369918513056329728 

^ ' ^ ~~ 765014898375 ' ^ ' ' ~ 79496721732575 ' 

_ 6471153009519628976971776 _ 2239521974640025214976 

^'' ) ~ 529050683130286625 ' 397483608662875 ' 

512 n 10 



W(n) 



81 (9 n 2 - 64) 2 (81 n 2 - 1024) (81 n 2 - 64) 4 (81 n 2 - 256) 4 



x -447359221655207230383849472 + 3163532546828444207717285888 n 2 



24 



-9812275706440214903775035392 n 4 + 17549420408786506395496218624 n 6 
-20019188843798683542274179072 n 8 + 15231566821196295158466871296 n 10 
-7874351861115834918827458560 n 12 + 2775326655204987593042165760 n 14 
-660250519084024545000996864 n ie + 103213969217030305334486016 n 18 
-10045634813323035262758792 n 20 + 545473369543317738837741 n 22 
-12459805959909788031573 n 24 ) , (A.3) 



_ 1670626488940208128 525333494026541203456 
^ ' ~ 485563688671875 ' ^ ' " 26299512173025 ' 

_ 6910180342493957842421202407424 
^ ^ ~ 2280028384216811071375 ' ^ ' ' 



B Octet potential 

In this appendix we briefly discuss the role played by the octet coulomb potential and 
the error done neglecting it. 

Frequently the octet potential is neglected in the octet propagator when the gluon corre- 
lation time is taken into account || § (see |7] where the octet potential is taken into account 
being solved the differential equations numerically). Therefore, we find pretty interesting to 
consider the case for a general octet potential and take the limit N c — > oo in our results. 
Thus, we recalculate H and W with | = Af2 1 _ 1 = r for a general r. That is, 

T , 1 & s 1 C F a s 

2N C r N£ — 1 r 

H(n,l) -> H oct (n,l,r) , 

W(n)^W oct (n,r). (B.l) 

We obtain for r = (N c — > oo) 

H oct (l, 0, 0) = 198 ~ 1.18#(1, 0) , H oct (2, 0, 0) = 312576 ~ 1.32#(2, 0) , 

# oct (2, 1, 0) = 183040 ~ 1.37#(2, 1) , H oct (3, 0, 0) = 22609206 ~ 1.38#(3, 0) , 
H oct (3, 1, 0) = 17058600 ~ 1.39#(3, 1) , H oct {3, 2, 0) = 7846956 ~ 1.39#(3, 2) , 
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, s 20548 . . , N 133338112 
W«*(1,0) = 1.19^(1), W oct (2,0) = 1.34^(2), 

«0) = -^~»(3). (B.2) 

We find they are always larger than the values obtained with N c = 3. We also see that 
the error is large though under control. Finally, we have also been able to obtain analytical 
expressions for W and H (with I = 0) for a general r. They read 



„, , x _ in /3 (1-n 2 ) 17 1-n 2 898 - 873n 2 - 25n 4 

W oct (n,r) = 288 n 10 (— ^ 1 + ^ / + ^ 5 ^ 

1 ' ^ V 2(l+r) 4 3(l + r) 3 72 (1 + r) 2 

17976 - 17741 n 2 - 235 n 4 ) n (120 - 274 n + 225 n 2 - 85 n 3 + 15 n 4 - n 5 



864 (1 + r) 17280 (-A + n + nr) 

-3 + n) 2 (-40 + 94n-77n 2 + 26n 3 -3n 4 ) 



1440 (-3 + n + nr) 2 



2520 - 9186 n + 13531 n 2 - 10240 n 3 + 4160 n 4 - 854 n 5 + 69 n 6 ) 



4320 (-3 + n + nr) 



-2 + n) 4 (-3 + 4n-n 2 ) (2 - n) 3 (108 - 75 n - 56 n 2 + 23 n 3 ) 



+ 



24 (-2 + n + nr) 4 432 (-2 + n + nr) 3 

-2 + nf (-360 - 4743 n + 8026 n 2 - 3317 n 3 + 394 n 4 ) 



4320 (-2 + n + nrf 
■10080 + 99624 n - 142806 n 2 + 42615 n 3 + 22325 n 4 - 13479 n 5 + 1801 n 6 ) 



17280 (-2+n + nr) 
2 (2-n) (-1 +n) 4 n (-1 + n) 3 n (126 - 53 n — 5 n ) 

3(-l+n + nr) 4 27(-l+n + nr) 3 

(-1 + nfn (8070 - 1103 n- 1612 n 2 + 73 n 3 ) 



864 (-1 + n + nrf 
n (-121380 + 87068 n + 87765 n 2 - 60855 n 3 + 9615 n 4 - 2213 n 5 ) 



8640 (-1 + n + nr) 
2 n (1 + nf (2 + n) n (1 + n) 3 (-126 - 53 n + 5 n 2 ) 



3(l+n + nr) 4 27(l + n + nr) 3 

n (1 + nf (-8070 - 1103 n + 1612 n 2 + 73 n 3 ) 



864(1 + n + nr) 2 
n (-121380 - 87068 n + 87765 n 2 + 60855 n 3 + 9615 n 4 + 2213 n 5 ) 



8640 (1 + n + nr) 
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(2 + n) 4 (3 + 4n + n 2 ) (2 + nf (-108 - 75 n + 56 n 2 + 23 n 3 ) 

24(2 + n + nr) 4 432 (2 + n + nr) 3 

(2 + n) 2 (-360 + 4743 n + 8026 n 2 + 3317 n 3 + 394 n 4 ) 

4320 (2 + n + nr) 2 
(10080 + 99624 n + 142806 n 2 + 42615 n 3 - 22325 n 4 - 13479 n 5 - 1801 n 6 ) 

17280 (2 + n + nr) 
(3 + n) 2 (40 + 94n + 77n 2 + 26n 3 + 3n 4 ) 

1440 (3 + n + nr) 2 
(2520 + 9186 n + 13531 n 2 + 10240 n 3 + 4160 n 4 + 854 n 5 + 69 n 6 ) 

4320 (3 + n + nr) 
n (120 + 274 n + 225 n 2 + 85 n 3 + 15 n 4 + n 5 ) s 
17280 (4 + n + nr) ' ' 



, , s/ 72(-l+n 2 ) 228 (-1 + n 2 ) 2 (-214 + 229 n 2 - 15 n 4 ) 
# orf (n,0,r) = n 8 (— \ + £ + 



(1 + r) 3 {l + ry 1 + r 

-1 + nf (-24 + 26 n - 9 n 2 + n 3 ) 2 (-2 + n) 3 (3 - 4 n + n 2 ) , . 

+ » T3 ( B - 4 ) 



2 (-3 + n + nr) (-2 + n + nr 

'-2 + nf (3 + 20 n - 31 n 2 + 8 n 3 ) (12 - 1 12 n + 117 n 2 + 18 n 3 - 45 n 4 + 10 n 5 



(-2 + n + nr) 2 -2 + n + nr 

32 (-2 + n) (-1 + n) 3 n 8 (-1 + n) 2 n (-22 + n + 5 n 2 ) 



+ 



(— 1 + n + n r) 3 (— 1 + n + nrf 

n (578 - 303 n - 539 n 2 + 255 n 3 + 9 n 4 ) 32 n (1 + n) 3 (2 + n) 



2 (-1+n + nr) (i +n + nr y 

8 n (1 + n) 2 (22 + n - 5 n 2 ) n (578 + 303 n - 539 n 2 - 255 n 3 + 9n 4 

(1 + n + nr) 2 2 (1 + n + nr) 

2 (2 + n) 3 (3 + 4n + n 2 ) (2 + n) 2 (3 - 20 n - 31 n 2 - 8 n 3 ) 



+ 



(2 + n + nr) 3 (2 + n + nr) 2 

■12 - 112 n - 117n 2 + 18 n 3 + 45 n 4 + 10 n 5 ) (1 + nf (24 + 26 n + 9 n 2 + n 3 



2 + n + nr 2 (3 + n + nr) 

( |B.3|) and ( p.4j) do not have a unique limit when n — > integer, r — > 0. The proper order is 



lim lim 

T — >Q n — 'integer 
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